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ABSTRACT 

Context. The elements technetium and lithium are two important indicators of internal nucleosynthesis and mixing in late type stars. 
Studying their occurrence and abundance can give deep insight into the structure and evolution in the late phases of the stellar life 
cycle. 

Aims. The aims of this paper are: 1) to revisit the Tc content of a sample of oxygen-rich asymptotic giant branch (AGB) variables 
and 2) to increase the number of such stars for which the Li abundance has been measured to provide constraints on the theoretical 
models of extra-mixing processes. 

Methods. To this end, we analysed high-resolution spectra of 18 sample stars for the presence of absorption lines of Tc and Li. 
The abundance of the latter was determined by comparing the observed spectra to hydrostatic MARCS model spectra. Bolometric 
magnitudes were established from near-IR photometry and pulsation periods. 

Results. We reclassify the star V441 Cyg as Tc-rich, and the unusual Mira star R Hya, as well as W Eri, as Tc-poor. The abundance 
of Li, or an upper limit to it, was determined for all of the sample stars. In all stars with Tc we also detected Li. Most of them have 
a Li content slightly below the solar photospheric value, except for V441 Cyg, which has ~1000 times the solar abundance. We also 
found that, similar to Tc, a lower luminosity limit seems to exist for the presence of Li. 

Conclusions. We conclude that the higher Li abundance found in the cooler and higher luminosity objects could stem from a 
production mechanism of Li operating on the thermally pulsing AGB. The stellar mass might have a crucial influence on this (extra 
mixing) production mechanism. It was speculated that the declining pulsation period of R Hya is caused by a recent thermal pulse 
(TP). While not detecting Tc does not rule out a TP, it indicates that the TPs are not strong enough to drive dredge-up in R Hya. 
V441 Cyg, on the other hand, could either be a low-mass, intrinsic S-star that produced its large amount of Li by extra-mixing 
processes, or an intermediate-mass star (M > 4Mq) undergoing Li production due to hot bottom burning. 

Key words, stars: AGB and post-AGB - stars: abundances - stars: evolution - stars: interiors - stars: individual (R Hya, V441 Cyg) 



1. Introduction 

The asymptotic giant branch (AGB) phase is the last stage of 
evolution with nuclear burning for stars of low and intermediate 
mass (1 - 8M0). In the most luminous part of the AGB, the be- 
haviour of an AGB star is characterised by the so-called thermal 
pulses (TP), thermal instabilities of the He shell accompanied 
by changes in luminosity, temperature, pulsation period, and in- 
ternal structure (see e.g. |Busso et al.| |1999| [Herwig 2005] for 
reviews). Between the repeated events of explosive He-burning, 
heavy elements can be produced via the "slow neutron capture 
process" (s-process, see e.g. |Wallerstein et al.| [1997) in the re- 
gion between the hydrogen- and the helium-burning shells. The 
processed material may then be brought to the stellar surface 
by the convective envelope that temporarily extends to these 
very deep layers. This mixing event is called the third dredge- 
up (3DUP) and is the cause of the eventual metamorphosis of an 
oxygen-rich M-star into a carbon-rich C-star. 



* Based on observations at the Very Large Telescope of the European 
Southern Observatory, Cerro Paranal/Chile under Programme 65.L- 
0317(A), and at the Calar Alto Observatory Centra Astronomico 
Hispano Aleman, Calar Alto, Spain 



The element technetium (Tc, Z = 43) is a well-established 
diagnostic tool in the study of the thermally pulsing AGB (TP- 
AGB) phase. Because of the short half life time of the longest 
lived Tc isotope produced in the s-process ("Tc, T\pi = 2.1 x 
10 5 yr), any Tc that is detected in the star's atmosphere must have 
been produced in-situ. This requirement has been used to study 
the evolution of AGB stars ever since the pioneering discovery of 
Tc absorption lines in the spectra of late type giants by Merrill 
( 1952b| l. Several classes of objects have been investigated for 
the presence of Tc in the past (e.g.|Dominy & Wallerstein| |l986; 
|Litt e et al.| |1987| |Van Eck & Jorissenj |1999| l. A key "question 
in this research is if and when in the evolution on the AGB, 
in terms of luminosity, the 3DUP sets in. This question was 
addressed by |Lebzelter & Hroii] ( |1999| |2003[ ) for a sample of 
solar-neighbourhood semi-regular pulsating AGB stars and for 
a larger sample of Galactic disk AGB stars, including several 
Mira variables. Recently, [Uttenthaler et al.| ( |2007a| l has investi- 
gated a sample of Galactic bulge AGB variables, for which the 
distance and hence luminosity is more precisely known than for 
their Galactic disk counterparts. The general conclusion of these 
works is that the lower luminosity limit for 3DUP to set in, as 
predicted by theoretical evolutionary models of the AGB, is con- 
sistent with the lowest luminosity Tc-rich stars. However, being 
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above this limit is not sufficient for a star to have Tc in its at- 
mosphere ( |Lebzelter & Hron| |2003| |Uttenthaler et al.| [2007a i. 
In the present paper we revisit the Tc content of a sample of 
oxygen-rich AGB stars, with important consequences for a few 
individual stars. 

The element lithium (Li, Z = 3) is very fragile and 
temperature-sensitive. Once it is mixed into stellar layers with 
temperatures exceeding 3 x 10 6 K, it is destroyed via proton cap- 
tures to form helium. Its fragility makes Li an important diag- 
nostic tool for stellar evolution (Rebo lo) |1991| l, and also has a 
high significance in the determination of cosmological parame- 
ters ( |Spite & Spite||1982||Korn et al.||2006) . The standard model 
(e.g. |Michaud & Charbonneau||1991| l predicts a strong reduction 
in the Li surface abundance after the first dredge up to values 
log e(Li) < 1 .5, or even less if depletion effects on the main se- 
quence are considered. Most giants for which a Li abundance has 



been determined (see a summary by Michaud & Charbonneau 
|1991) > agree with this result. 

Despite its fragility, Li can be produced rather than destroyed 
under certain condition in stars via the so-called Cameron- 
Fowler or 7 Be-transport mechanism ( 3 He(a,y) 7 Be(<?~, v) 7 Li; 
|Cameron & Fowler] |1 97 1) . The production of Li is always con- 
nected to mixing mechanisms at or below the bottom of the con- 
vective envelope, and requires 3 He to be present in the stellar 
atmosphere. One stage in stellar evolution where Li production 
occurs is the luminosity bump of the red giant branch (RGB) 
when the H-burning shell erases the molecular weight disconti- 
nuity left behind by the first dredge up. A small fraction of giants 
in this phase are known to be very Li-rich ( de la Reza et al. [1997 ; 

2000| l. The mixing process below 



Charbonnel & Balachandran 



the convective envelope has been named extra-mixing or cool 



bottom processing (CBP; Wasserburg et al. 1995 i. Its effects 



are also reflected in the carbon isotopic ratio measured from the 
spectra of giant stars ( |Gilroy & Brown|[l991| l, as well as in other 
isotopic ratios measured in meteoritic pre-solar grains of RGB or 
AGB origin ( |Nollett et al.|[2003] l. 

The physical cause of the extra-mixing below the convec- 
tive envelope is hitherto unknown. Eggleton et al. ( 2006| l, on 
the basis of a 3D simulation, find that Rayleigh-Taylor insta- 
bilities below the convective envelope can develop from the in- 
version of the mean molecular weight gradient induced by 3 He 
burning. Alternatively, Charbonnel & Zahn ( 2007 ) suggest that 
the double-diffusive mechanism called "thermohaline instabil- 
ity" should be at play. In principle, both these phenomena might 
also occur on the AGB, though detailed models are just emerging 
( |Cantiello & Langerj [2008) 1. Finally, |Bus"so et al.| ( |2007] l explore 
the possibility that circulation of partially processed matter can 
be accounted for by the magnetic buoyancy induced by a stellar 
dynamo operating on the RGB and on the AGB. That CBP is 
probably active on the AGB has been recently demonstrated by 
lUttenthaler etal~| ( 12007b| l and |LebzelterlM^ ( |20()8] i. 

Among the AGB stars, those with carbon-rich chemistry (C- 
type, C/0>1 due to the dredge-up of freshly synthesised 12 C) 
have been studied more intensively for their Li content than 
those with oxygen-rich chemistry (M-type, C/0<1) in the past. 
This may partly have a historical reason because of the early de- 



tection of carbon stars with a strong Li resonance line (McKellar 
1940[ |Torres-Peimbert et al.| |1964| l, but also an observationa 
reason. Although the molecular line blanketing mainly caused 
by the CN and C2 molecules in the spectral vicinity of the Li 
line in C-stars is strong, that of the TiO molecule in oxygen-rich 
stars of the same effective temperature is even more severe. The 
lower limit of Li abundance that can be measured is thus higher 
for O-rich stars than for C-rich stars. In both cases, a Li abun- 



dance determination can only be achieved with the help of spec- 
tral synthesis techniques including large numbers of molecular 
lines. Thus, the counterparts to the catalogues of Li abundance 
measurements in C-stars ( |Denn et aL] 1 1 99 1 [ |Boffin et al.[[ l993) 
include only relatively small tables of M-type giants ( |Luck & 
Lambert, 1982). Even in this last work, a number of the observed 
stars are supergiants and non-variable giants that are rather in the 
RGB than in the AGB phase of evolution. 

There are, however, M-type AGB stars that are known to 
be very rich in Li. This type of star has been observed in the 
Magellanic Clouds ( Smith et al. 1995|l, and very recently also 
in the Milky Way ga" ' 



axy ( Garcfa- 



lernandez et al. 11200711. These 



objects are intermediate-mass (M > 4 - 8Mq) TP- AGB stars 
with a convective envelope penetrating the H-burning shell, thus 
nuclear-burning occurs partly under convective conditions. This 



process is generally known as hot bottom burning (HBB; Iben 
|1973| ). Under these conditions, carbon, on the one hand, is con- 
verted into nitrogen, keeping the C/O ratio below 1 , and on the 
other hand 7 Li is produced very efficiently by the Cameron- 
Fowler mechanism. The Li abundance on the surface of HBB 
stars can be higher than the initial abundance of the interstel- 
lar medium from which they formed. We do not intend to deal 
with this kind of stars in the present paper, but rather with the 
less massive AGB stars that have not (yet) dredged up enough 
carbon to have C/0>1 and are not massive enough to remain O- 
rich due to the operation of HBB. Nevertheless, by inspecting 
their bolometric magnitudes, we will check whether our Li-rich 
objects belong to this group of stars. 

The correlation of Tc and Li was studied for a sample 
of Galactic S-type (C/O ^l) stars by |Vanture et~aT] ( |2007| i. 
|Uttenthaler et al.| ( |2007b[ ) present Li abundances and their cor- 
relation with Tc in a homogeneous sample of O-rich AGB stars 
located in the Galactic bulge. Here, we want to increase the sam- 
ple of O-rich AGB stars whose Li abundance has been measured, 
and study its correlation with the presence of Tc in this type of 
stars. These measurements will be useful for understanding the 
interaction between CBP mechanisms and 3DUP in AGB stars. 

The paper is structured in the following way. In Sect. [2] the 
sample and the observations are described and bolometric mag- 
nitudes of the sample stars derived; in Sect. [3] we present the 
analysis of the spectra, results concerning the presence of Tc and 
the abundance of Li; in Sect. [4] the results are discussed in view 
of a possible enrichment mechanism of Li, and we take a closer 
look at two of the sample stars; Sect. [5] summarises the work. 



2. Sample and observations 

The sample used in this study originally goes back to the search 
for Tc in Galactic -disk AGB variables in the solar neighbour- 
hood, based on parallaxes measured with the Hipparcos satellite 
( |Lebzelter & Hron |2003[ ). Primary targets for that study were 
those stars that had not been searched for Tc in the past and 
objects without a definite decision on their Tc content. The orig- 
inal spectra of that study are reanalysed in this work. An excep- 
tion here is the star V441 Cyg, whose observation and analy- 



sis was initiated by the work of Vanture et al. (2007). A refined 



method for detecting Tc described in Sect. |3.1| has led to some 
re-classifications in the sample, which we describe in the present 
study. 

Table[T]lists the 18 objects studied in this project with a num- 
ber of important characteristics. The spectral type is a kind of 
"mean" value for what is found in the literature (using Vizier; 
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Ochsenbein et al. 2000 1^ The carbon star TW Hor is an excep- 
tion in this sample. We include it because it was observed in the 
same UVES run as most of the other sample stars. According 
to the spectral type designations, the stars W Cet and V441 Cyg 
have a notable enhancement of s-elements in their atmospheres. 
V441 Cyg is described in the investigation of Galactic S-stars 
by |Vanture et al.| p007 ), who identify this star as the only one 
in their sample showing Li in its spectrum, but not Tc. We come 
back to this star in Sect. 4. 1 where we challenge this classifica- 
tion. 

The variability type in column 3 of Table [T] is taken from 
the Combined General Catalogue o f Variable Stars (CGCVS), 
version no. 4.2 ( |Samus et al.| |2004| l. Also the pulsation period 
in column 4 is taken from the CGCVS, except for Y Scl, for 
which the period is taken from the All Sky Automated Survey 
catalogue (Pojmanski et al. 2005 | l, as well as RS Vir and S Vir, 
for which the period is taken from Temple ton et al.| ( [2005| l. A 
further exception is V441 Cyg. Because of the importance of this 
star for the further discussion of the results (see below), we made 
a check of the period listed in the General Catalogue of Variable 
Stars (375 d). We did a Fourier analysis of the last 3000 days 
of the AAVSO archive data of this st ar (|Henden| |2009] > using 
the program Period04 ( |Lenz & Bre ger 2005} . We found a quite 
well-expressed peak at 288 d and a possible second period at 
153 d. A typical visual light amplitude of 2™0 is found. 

The near-IR photometry listed in columns 5 and 6 of Table[T] 
is taken from the references listed in column 7. The derivation 
of the bolometric magnitude in column 8 is described in more 
detail in Sect. 12. II The IRAS colour listed in column 9 is defined 
as[12]-[25] = -2.51og(Fi 2 /F 25 ). 

Most of the targets listed in Table [T] were observed in 
July 2000 with the UVES spectrograph at ESO's Very Large 
Telescope under Programme 65. L-03 17(A). The resolution of 
these spectra is R — A/AA 50000, and the wavelength cov- 
erage is approximately 380 - 490 nm in the "blue arm", and 667 
- 1000 nm in the "red arm" of UVES. These observations are de- 

(2007a), and the spectra are analysed 



scribed in Uttenthaler et al. 



for Tc in Lebzelter & Hron 



(2003 ). We thus refer to these papers 



for more details. The only exceptions to this are the stars W Cet, 
which is not included in |Lebzelter & Hron| ( |2003| l, and V441 
Cyg, which was observed for the present study in service mode 
with the FOCES spectrograph (Pfeiffe r et al.||1998| ) mounted on 
the 2.2 m telescope of the Calar Alto Observatory in southern 
Spain on 23 June 2008. The observations were carried out using 
the CCD chip with 15 fim pixel size and a slit width of 100/im, 
which subtends 1'.' 15 on the sky. The chosen setup yields a spec- 
tral resolving power of R 65 000. The exposure time was two 
times 1800 s. Because of the high dynamic range of FOCES and 
the resulting reduction problems with the "red", "green", and 
"blue" flats, the spectra were reduced with a custom-made IDL 
routine. The signal-to-noise ratio (SNR) of the spectra is esti- 
mated to be in the range 6 - 10 at the wavelength of the Tc lines, 
and >100 at the wavelength of the 670.8 nm Li lin^j As shown 
by |Uttenthaler et al.|p007a} , an SNR of 6 - 10 is sufficient to 
decide on the presence of Tc based on a flux-ratio method (see 
next section), as long as the spectral resolution is high enough. 

The general appearance of the spectra is dominated by strong 
bands of the TiO molecule (except of the carbon star TW Hor, 
of course), clearly verifying their oxygen-rich nature. Some of 
the stars also show significant bands of the ZrO molecule, in- 
dicating that some s-process enrichment has taken place. Only 

1 http://webviz.u-strasbg.fr/viz-bin/VizieR 

2 All wavelengths in this paper are expressed in nm in air. 



stars that are identified as Tc-rich (see below) also show bands 
of ZrO, and conversely no Tc-poor star shows any clear sign of 
ZrO absorption bands. Thus, none of our sample stars is extrinsi- 
cally enriched in s-process elements by, e.g., mass transfer from 
a binary companion. However, the absorption bands of the LaO 
molecule cannot be clearly identified in any of our sample stars, 
as is sometimes the case in S-type stars. 

2.1. Bolometric magnitude determination 

To derive absolute /^-magnitudes, we used the logP-K relations 
from lta et al.| ( |199"4"| ). For stars classified as Miras we used Ita's 
relation C including their separation for targets with J - K < 1 .4 
and J - K > 1 .4. For semi-regular variables we used sequence 
B or C. Naturally, this approach includes some uncertainties. 
In their Table 3, |Ita et al.| ([1994) list a typical width of about 
07*2 of these sequences. This, however, is much less than the 
uncertainties of absolute Zf-magnitudes derived from parallax 
measurements, which are only available for a small number of 
stars in our sample. A larger error comes from the uncertainty of 
attributing a star to a given pulsation sequence. Several authors 



agree ( |Wood etaL] |1999 Kiss & Beddi ng 2003 Lebzelter et al. 
2005) that large amplitude variables, classically dubed Miras 
are found on sequence C, which corresponds to fundamental 
mode pulsation. The risk of a misclassification will thus be low 
for the Miras. It is much more difficult to attribute the inhomoge- 
neous group of the five semi-regular sample stars correctly. They 
can be located on any sequence from A to C, and the definition 
of sequences A, B, and C is less clear. As selection criterion, 
we chose the light amplitude following the extensive analysis of 
Kiss & Bedding (2003) for the LMC. As on average the ampli- 



tude becomes smaller from sequence C to A, we attributed the 
non-Miras V441 Cyg, Y Scl, and EV Vir to sequence C, TW 
Hor and ER Vir to sequence B. It is clear that this classification 
and the derived luminosities have to be taken with some caution. 
However, because of the small separation of sequences B and C 
compared to the difference to sequence C, any misclassification 
will not change the overall result substantially. 

To derive absolute bolometric magnitudes from the absolute 
^-magnitudes, we need to apply a bolometric correction. We 
adopt the approach outlined in Mekul (2005) using two different 
relations between blackbody fits and J — K depending on the 
K - L colour. The inclusion of a second near infrared colour 
allows for much better treatment of stars at the cool end of the 
AGB than other formulas found in the literature. Details will be 
given elsewhere (Kerschbaum et al., in preparation). 

We used near-infrared colours from the literature (see 
Table [T). The J and K band magnitudes were taken from 
I Catchpole et~aT[ ([19791), |Fouquet et al] ( [1992] ), [Kerschbaum &| 
|Hron| ( |1994} and Kerschbaum| ( 1995 ), respectively. All infrared 
measurements were converted to the 2MASS system using the 
relations given in Carpenter (2001 ). Where more than one mea- 
surement was found we used an average value. While the K- 
magnitudes used to derive the absolute bolometric magnitudes 
were calculated from a star's period, the bolometric correction 
is based on individual near infrared measurements and thus af- 
fected by stellar variability. The extensive material presented by 
Catchpole et al.| JT979 ) allow estimation of the typical uncer- 
tainty of the bolometric correction from the variations in the 
J - K colour. For Miras, this uncertainty can be up to 0?5. In 
a few cases no K—L colour was available. We then used an aver- 
age of the two relations. In all cases the error introduced by this 
is less than O^l. Interstellar extinction has not been taken into 
account for these rather nearby objects. 



4 



Uttenthaler & Lebzelter: Correlation between Tc and Li in AGB variables 



Table 1. I mportant characteristics of th e s ample stars. Referenc es for near infrared data: 1 : Catchpole et al. ( 1979 1; 2: Fouquet et al. 
( fT992l ; 3: |Kerschbaum & Hron| (19941 ); 4: |Kerschbaum| | fl995) . 



Name 



Sp. type Var. type P [days] K [mag] J - K [mag] Ref 



[12] -[25] 



UCet 


M4 


M 


234.76 


2?73 


1?11 


1 


-4™60 


-1.14 


WCet 


S7 


M 


351.31 


2? 14 


1?23 


1 


-5?07 


-1.32 


V441 Cyg 


M4S 


SRA 


288.00 


1"?69 


T39 


3 


-5?66 


-1.28 


W Eri 


M7 


M 


376.63 


r?78 


1?42 


1 


-5^08 


-0.90 


S Gru 


M7 


M 


401.51 


o: , 54 


1?22 


1 


-5P35 


-0.87 


RHor 


M7 


M 


407.60 


-CT88 


1?24 


1 


-5?35 


-0.92 


THor 


M5 


M 


217.60 


3™34 


1™13 


1 


-4?44 


-1.15 


TW Hor 


C5(N) 


SRB 


158.00 


a: 1 19 


l'M6 


3 


-4?89 


-1.03 


RHya 


M7 


M 


388.87 


-2™52 


1?24 


1 


-3?28 


-1.08 


RUHya 


M6 


M 


331.50 


l m 60 


I'M 8 


1 


-5?07 


-0.74 


YLib 


M5 


M 


275.70 


3™16 


1?21 


1 


-4™72 


-0.73 


UMic 


M6 


M 


334.29 


l m 80 


1™22 


1 


-5™00 


-0.69 


YScl 


M5 


SRB 


79.36 


0P35 


1™27 


2 


-4?24 


-0.94 


TTuc 


M4 


M 


250.30 


2 m 96 


1?18 


2 


-4™58 


-1.08 


ER Vir 


M4 


SRB 


55.00 


l'M4 


no 


4 


-4™03 


-1.41 


EV Vir 


M4 


SRB 


120.00 


1?55 


ri5 


4 


-5™08 


-1.34 


RS Vir 


M7 


M 


353.35 


ri5 


1™31 


1 


-5^09 


-0.56 


S Vir 


M7 


M 


378.01 


0?30 


1?27 


1 


-5P21 


-0.98 



-5.4 



8-5.2 



o 

T5 



-4.8 



o 
n 

2-4.6 



-4.4 



-4.4 



-4.6 



-4.8 -5 
Mbol (this work) 



-5.2 



-5.4 



Fig. 1. Comparison of M\, Q \ values derived by out method com- 
pared with values using the formula given in Guandalini & 
Busso (2008). The full line indicates the 1:1 relation. Good 
agreement can be seen. 



Recently, |Guandalini & B usso ( 2008 ) have published a new 
set of bolometric magnitudes for mass-losing AGB stars. As a re- 
sult they give a log P-M\, \ relation for Miras (see their Fig. A. 1). 
Using this relation we calculated Mt, i for the Miras in our sam- 



ple. In Fig.[Tlwe compare the results derived from Guandalini & 



Busso (2008 ) with our values of the absolute bolometric magni- 
tude for the stars in common, showing good agreement between 
the two approaches. 

3. Analysis and results 

3.1. Analysis with respect to Tc 

The radial velocity (rv) is the first quantity that we determined 
from our spectra. It was measured from the TiO band head at 
705.4 nm using a cross-correlation technique, with a synthetic 
MARCS spectrum as template. Sect. 3.2.1 gives more details on 
the used MARCS models. We chose this TiO feature because 
the SNR of the observed spectra is much higher at this wave- 



length than at the wavelength of the Tc lines, and the TiO fea- 
tures are formed in the stellar atmosphere close to the atomic 
(Tc) lines. For TW Hor, the rv was measured from atomic lines 
in the vicinity of the Tc lines, due to the absence of TiO bands 
in carbon stars. The precision of the rv measurement is of the 
order of lkms 1 or better. The differences to literature values 
are usually small (a few kms~'), and can be fully explained by 
rv variations caused by the stellar pulsation. 

We corrected the observed spectra by the measured rv shift 
to check them for the Tc lines using the flux-ratio method intro- 
duced by |Uttenthaler et aL ( 2007a| >. In this method, the mean flux 
in a small wavelength interval around a quasi-continuum point 
close to the Tc line is divided by the mean flux in a small wave- 
length interval centred on the Tc line itself. This ratio can be cal- 
culated for all of the seven identified Tc lines (398.497, 403.163, 
404.911 and 409.567, 423.819, 426.227, and 429.706 nm), and 
plots of the ratios from different lines clearly separate the Tc- 
poor from the Tc-rich stars. If the Tc lines are absent, the ratio 
will be close to 1.0, but significantly more than 1.0 if they are 
present. Figure [2] shows such a plot for the lines at 423.819 and 
426.227 nm for our sample stars; plots involving other Tc lines 
are very similar and reveal the same Tc-rich stars. The Tc-rich 
stars that clearly separate from the Tc-poor ones are identified in 
the plot with their name. The error bars on the symbols in Fig. [2] 
were directly derived from the SNR of the observed spectra. 

We want to briefly explain the advantages of the flux-ratio 
method over previous methods of deciding on the presence of 
Tc. The method previously applied uses the central wavelength 
of the Tc line blend at 429.706 nm (e.g. |Little et al.| [1987] ; [v5n| 
|Eck & Jorissen| |1999| and references therein). This Tc line is 
blended with a fine of chromium. If Tc is present, the central 
wavelength is shifted to a somewhat longer wavelength than in 
the absence of Tc. However, the shift is not large, thus the sepa- 
ration between Tc-poor and Tc-rich stars is not very large either. 
Every uncertainty of the local wavelength scale has a direct im- 
pact on the central wavelength of the blend, thus adding uncer- 
tainty to the classification with respect to Tc. Due to the stellar 
pulsation and the velocity fields in the atmosphere, lines formed 
at different layers have a slightly different rv shift. This shift 
also depends on the pulsation phase at which the observations 
were done. The result is that intermediate cases are classified 
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3.5 



3.0 



ffi 2.5 



2.0 



£ 1.5 



1.0 



0.5 




V441 C yg * 

S Gru 



1.0 



1.2 



1.4 1.6 1. 

cont. / line Tc423.8 



2.0 



2.2 



Fig. 2. Ratio of the continuum flux to line core flux for the two 
Tc lines at 423.8 19 and 426.227 nm. The quasi-continuum points 
were taken in the range 423.940 - 423.970nm and 426.130 - 
426. 150 nm, respectively, the line core flux was measured in the 
range 423.805 - 423.835 nm and 426.220 - 426.240 nm, respec- 
tively. While the Tc-poor stars gather around ratios 1.0 (lower 
left corner), the Tc-rich stars are displaced to the upper right by 
the Tc line absorption. The dashed lines represent the upper lim- 
its of the ratios found for the Tc-poor bulge stars analysed in 
lUttenthaler eTaT1 ( |2007a| l. 



as "Tc doubtful", "Tc possible", and "Tc probable". Figure [2] 
shows that, with our flux-ratio method, the gap between Tc-poor 
and Tc-rich stars is so large compared to the uncertainty that no 
room is left for these intermediate cases. The uncertainty of the 
wavelength scale has only a very small impact on this method. 
We can thus very clearly decide on "Tc yes" or "Tc no". We 
conclude that the flux-ratio method delivers more reliable results 
than the central wavelength method. 

Table [2] summarises the radial velocity shifts measured from 
the spectra and the updated Tc classification. With respect to pre- 
vious classifications, there are three changes: V441 Cyg is Tc- 
rich (|Vanture et ak) |2007| Tc-poor), W Eri is Tc-poor fLebzelter| 
|& Hron|"[2 003 l Tc possible), and R Hya is Tc-poor as well 
( [Lebzelter & Hron) [20031 Tc-rich). 



3.2. Analysis with respect to Li 

The abundance of Li in stars is usually measured from the strong 
resonance doublet of neutral Li at 670.8 nm. The line strength 
is quite sensitive to the effective temperature of the star, thus a 
good knowledge of this stellar parameter is required for reliable 
abundance determination. In the following we describe in some 
detail how the temperature determination was accomplished for 
the present sample of stars. 



3.2.1. Effective temperature determination 

For determining the effective temperatures of the sample stars, 
we applied a x 2 minimisation method to band heads of the TiO 
molecule using a small grid of synthetic model spectra. In the 
optical range, TiO bands are the dominant spectral features of 
cool oxygen-rich giants, and they are very sensitive to the stellar 



temperature. We thus followed a procedure similar to Garcia 



Table 2. Quantities measured from the spectra: radial velocity, 
Tc content, effective temperature, and Li abundance of the sam- 
ple stars. 



Name 


rv [km s '] 


Tc? 


7eff 


loge(Li) 


UCet 


-19.0 


no 


3300 


< 0.0 


WCet 


+5.8 


yes 


3100 


+0.8 


V441 Cyg 


-61.7 


yes 


3400 


+4.6 


WEri 


+17.8 


no 


2950 


+0.4 


SGru 


-18.4 


yes 


2900 


+0.8 


RHor 


+56.0 


yes 


2900 


+0.6 


THor 


+44.9 


no 


3300 


< 0.0 


TWHor 


+ 14.5 


yes 


3000 


-1.5 


RHya 


-15.7 


no 


3100 


<0.0 


RUHya 


+0.4 


no 


3300 


< 0.0 


YLib 


+9.3 


no 


3100 


+0.9 


UMic 


-61.0 


no 


2900 


+0.6 


YScl 


+30.0 


no 


3200 


+ 1.0 


TTuc 


-37.9 


no 


3300 


<0.0 


ER Vir 


+8.5 


no 


3500 


< 0.0 


EV Vir 


-44.5 


no 


3500 


< 0.0 


RS Vir 


-15.9 


no 


2900 


+ 1.0 


S Vir 


+7.2 


yes 


2800 


+0.5 



abundances in massive Galactic AGB (OH/IR) stars. We used the 
same modified^ 2 test as these authors to fit model data Y synt ^i to 
observed data Y i, s j, which is expressed in the form 



X 



S 

i=1 



[Y ob s,i ~ i{X\ . . . Xm)\ 



(i) 



Hernandez et al. ( 2007) 1, who apply this method to measure Li 



with N the number of data points and M the number of free 
parameters. The synthetic spectra are based on COMARCS at- 
mospheric models. COMARCS is a revised version of MARCS 
( J0rgensen et al. |1992| l with spherical radiative transfer routines 
from [Nordlund ( 1984]> a nd new opacity data from the COMA 
program ( Aringer|2000| i. The only model parameters that we 
varied in this procedure were the effective temperature and the 
metallicity. The other stellar parameters have a very limited im- 
pact on the spectral appearance, hence on the results of this 
analysis (see also Garcfa-Hernandez et al. 2007 for a more de- 
tailed discussion). For the surface gravity logg [cms 1 ], micro- 
turbulent velocity and the mass M, we chose values typical 
for red giant stars of 0.0, 3 kms -1 , and 1 Mq. Since the precise 
C/O ratio cannot be determined from our optical data (near-IR 
spectra would be required to do so), we fixed it at the solar value 
of 0.48, although in principle it might deviate from that value, 
particularly for the Tc-rich sample stars. Once the C/O ratio ap- 
proaches 1, the chemistry and thus the dominant features in the 
spectrum will change significantly. Two stars of our sample are 
classified as MS or S-type (W Cet and V441Cyg), and only for 
these objects do we expect that the assumption of a constant C/O 
might deteriorate the results to some degree. We come back to 
this point in the next section. 

We chose three metallicity values for the grid of model atmo- 
spheres: 0.3, 1.0, and 1.6 times the solar metallicity. The temper- 
ature ranges between 2600 and 3600 K in steps of 100K, with 
intermediate steps at 2950 and 3150K for the solar metallicity 
models. From this grid of model atmospheres, synthetic spectra 
were calculated in the two spectral pieces 668 - 674 nm (cov- 
ering the Li line at 670.8 nm and the TiO y(l,0)Rc and TiO 
y(l,0)Rb band-heads) and 700 - 710nm (covering the strong 
TiO y(0,0)Ra band-head at ~ 705 nm and the TiO y(0,0)Rb 
band head). Atomic lines from the VALD database ( Kupka et al. 
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1999) and TiO lines from the list o f|Schwenkel ( |1998| ) were taken 
into account in the spectral synthesis. Other molecular line data 
of less impact on the spectral ranges under consideration were 
also taken into account; these data are summarised in Table 1 of 
ICristallo et alT] ( |2007l >. 

Prior to the x L minimisation fitting, the model spectra were 
convolved with a Gaussian to a resolution of R = A/AA = 50 000 
(except for V441 Cyg, for which the spectra are convolved to 
R = 65 000). An additional macro-turbulence of 3kms _1 was 
assumed. This proves an acceptable value for all except two 
stars, namely V441 Cyg and Y Scl, for which a higher value 
of 6kms~ 1 had to be adopted to make a satisfactory fit to the 
spectra. 

Although the metallicities adopted for the atmospheric 
model grid cover a quite wide range, it turned out that the tem- 
perature at which y 2 becomes minimal does not strongly depend 
on the metallicity. The x 2 minimum for the spectral piece in- 
cluding the Li line is quite flat for many stars, thus the adopted 
temperature was mostly fixed by the second piece between 700 
- 710 nm. For the spectral synthesis to determine Li abundances 
(see next section), we therefore adopted the model with solar 
metallicity, with a temperature that satisfactorily fits both spec- 
tral pieces. There is some indication in the literature that solar 
metallicity is a good assumption for solar neighbourhood AGB 



stars (e.g. Vanture & Wallerstein 



2007) l. After finishing this procedure, we found that the adopted 



2002 Garcfa-Hernandez et al. 



effective temperature correlates well with the (J - K) colour, as 
can be expected: cooler stars have higher (redder) (J - K) colour 
than hotter stars. The adopted temperatures are summarised in 
the fourth column of Table [2] We estimate that the precision of 
the temperature determination is ±100K. 

3.2.2. Li abundances 

With this temperature, spectral synthesis calculations of the Li 
line doublet were performed to measure the Li abundance, us- 
ing the MARCS models with solar metallicity. The log gf val- 
ues from the VALD database were adopted for the 670.8 nm 
Li doublet (log gf = -0.009 and -0.309 for the component at 
the shorter and longer wavelength, respectively). The Li abun- 
dance was varied in steps of 0. 1 dex, and the observed-minus- 
calculated flux at the wavelength of the Li line was monitored 
to determine the Li abundance in the atmosphere of each sam- 
ple star. The abundance that left practically a zero flux differ- 
ence was finally adopted. The last column of Table [2] gives 
the measured Li abundances on the log e scale (log e(Li) = 
logN(Li)/N(H) + 12). The solar photospheric Li abundance on 
this scale is loge(Li) = +1.1. We also used our model spectra 
to verify that the observed absorption around the Li line cannot 
come from a line of cerium, which is observed to be a "sub- 
stitute" for Li in the hotter barium stars and post-AGB stars 
( |Lambert et aT|[T993||Reyniers et al.||2002] l. 

From the uncertainty of the temperature, the stellar param- 
eter with the strongest influence on the Li line strength, we es- 
timate an abundance uncertainty of ±0.4 dex. Taking the uncer- 
tainties of the other stellar parameters into account, which are 
harder to quantify, the total uncertainty on log e(Li) is estimated 
to be up to ±0.6 dex, being larger for the cooler stars. The rel- 
ative uncertainty within the sample, however, is probably less 
than this. 

A number of stars in the sample has no significant Li line 
absorption, their spectra are satisfyingly fit with negligible Li 
abundance in the spectral synthesis. It is difficult to state a pre- 
cise upper limit to their atmospheric Li abundance. In the vicin- 



ity of the Li line, positive and negative residuals are present in 
an observed-minus-calculated flux plot, due to imperfect spec- 
tral modelling. Besides errors on the adopted stellar parameters, 
neglecting velocity fields in the atmosphere (due to stellar pulsa- 
tions) by using hydrostatic models leads to persistent observed- 
minus-calculated flux residuals. Also, the line lists (mostly TiO) 
are certainly not without flaws at the high resolution of our spec- 
tra. The stars with no detectable Li line absorption are listed 
with log e(Li) < 0.0 in Table [2] an abundance that would cer- 
tainly not be detectable with our methods anymore. The carbon- 
star TW Hor was analysed by Kipper (2004), who determine a 
loge(Li) = -1.5 from another UVES spectrum of this object. 
We adopt this value of the Li abundance for TW Hor in Table [2 
We also compared our spectrum of TW Hor with the one usee 
by Kipper (2004) and found negligible variation between them. 
The low Li abundance of the carbon-star TW Hor would not be 
detectable in oxygen-rich stars of the same effective tempera- 
ture, because of the increased molecular absorption in M-type 
stars in this spectral region. In the following we focus on the 
oxygen-rich stars. Figure[3]shows sample spectra of two stars (R 
Hya and RS Vir) with very different Li abundances, along with 
a best-fitting model spectrum. The Li content of R Hya is below 
the detection threshold, whereas RS Vir is found to have Li at 
the solar photospheric level. 

The very high Li abundance of V441 Cyg is remarkable. 
As noted above, this star deviates from the trend towards higher 
(redder) J - K colours at lower temperatures: despite its rather 
high temperature estimate of 3400 K, a / — K = 1.39 is found. 
|Vanture & Wallerstein| pOO^ estimate T ef[ = 3100K from its 
spectral type. This low temperature is hardly compatible with 
our fit to the 700 - 710 nm spectral region. While there is some 
disagreement on the near-IR photometry (Wang & Chen 2002), 
the combination of red J-K but weak TiO bands could be caused 
by a high C/O ratio. From our models we found that a synthetic 
spectrum with C/O = 0.9 and a temperature of 3200 K also rea- 
sonably fits the 700 - 710nm region. We thus repeated the Li 
abundance measurement adopting C/O = 0.9 and T e ff = 3200 K 
and derived loge(Li) = 4.1, which is still a factor of 1000 



higher than the solar photospheric abundance! However, |Plez 
|et al.| ( |2003| l demonstrate that atomic and some molecular lines 
(such as from ZrO and LaO) strengthen even more when the 
C/O ratio approaches 1. To find the Li abundance that may be 
regarded as a lower limit, we constructed a MARCS model with 
a C/O ratio of 0.99 and T eff = 3400 K, but otherwise unaltered 
stellar parameters. Using this model, we find log e(Li) = 3.1 for 
V441 Cyg. The TiO bands are still well fit by the model at this 
high a C/O ratio. 

The equivalent width of the Li I 670.8 nm line in V441 
Cyg is approximately 730 mA, surpassing that of many of the 
super-Li rich stars identified in the Magellanic Clouds (Smith 
et al. 1995[ ), and the line core flux is ~ 6% of the local quasi- 



continuum level. Also, the Li line at 8 12.634 nm was identified 
in V441 Cyg. A spectrum synthesis of this line yielded a Li abun- 
dance lower by 0.8 dex (at solar C/O) than was found from the 
670.8 nm Li line. However, the spectral fit in the vicinity of the 
812.6 nm line was not as good as that around the 670.8 nm line, 
which caused problems in defining the quasi-continuum level. 

We repeated the Li abundance measurement with the C/O ra- 
tio increased to 0.99 for W Cet, the other S star in our sample. 
This exercise yielded log e(Li) = -2.0, which may again be re- 
garded as a lower limit for the Li abundance. The Li detection 
would still be significant, because the detection threshold will 
also decrease markedly with C/O approaching 1. The observed 
spectrum of W Cet is still incompatible with negligible Li con- 
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Fig. 3. Examples of observed spectra (black graph with squares) of two stars with very different Li abundances, along with a best- 
fitting model spectrum (grey graph). The dashed lines indicate the laboratory wavelengths of the Li doublet. 



tent, if C/O = 0.99 is adopted for the model atmosphere. For 
consistency, we keep the Li abundances as determined with the 
C/O = 0.48 models in Tableland in the analysis. 

We would also like to note that we did not identify significant 
LaO bands in any of our oxygen-rich stars, thus we may assume 
that indeed none of them has a C/O ratio above 0.99. 

Besides V441 Cyg, only one of our sample stars has pre- 
viously been investigated for its Li content, namely RS Vir 
( jGarcia-He rnandez et al.| |2007|l. In con trast to our measure- 
ment (Table |2 1, |Garcia-Hernandez et al. ( |2007| l find no signif- 
icant Li line absorption in that star, and state an upper limit of 
log e(Li) < 0.5 instead. Also, the temperature determinations of 
this star differ somewhat, 2900 K in the present work compared 
to 2700 K in |Garria-Hernandez et al.| ( |2007[ ). We identify two 
possible explanations for this discrepancy. 

First, two different TiO line lists were applied to derive stel- 
lar temperatures. One of them might yield systematically lower 
temperatures, which would decrease the measured Li abundance 
considerably. To investigate this possibility, we calculated model 
spectra using the TiO li ne list of |Plez (Q 998 ), as used by Garcia - 



( 2007 1, and repeated the^ 2 minimisation proce- 



Hernandez et al. 

dure (Sect. |3.2.1| >. Only the sequence of models with solar metal- 
licities was used for this. From this we find that the temperatures 
found using the two TiO line lists hardly differ from each other, 
with the Plez list delivering temperatures at most 50 K higher star certainly cannot be described perfectly with a hydrostatic 



than the Schwenke list. We conclude that the different line lists 
used cannot explain the discrepancy found in the temperature 
and Li content of RS Vir. From this experiment with the two 
TiO line lists, we also find that the list of |Schwenke| ([1998) 
yields slightly lower residuals than that of Plez|(1998i, partic- 
ularly in the vicinity of the 670.8 nm Li line. This might also 
explain why our detection limit is somewhat lower than that of 
IGarcia-Hernandez et al.| | |2007) . 

Second, the two spectra have been taken at different epochs, 
with different instruments. The temperature variation over the 
cycle of stellar pulsation could be responsible for the difference 
of 200 K found by the two studies. According to the AAVSO 
visual light curve, our spectrum was taken just 30 days before 
maximum, when the light curve was still sharply rising. Garcfa- 
|Hernandez et al.| ( [2007| l observed RS Vir during their fourth run 
between 21 - 24 February 1997. This date is just a few tens of 
days before a visual light maximum of the AAVSO light curve, 
too, but three pulsation cycles before our UVES observations. 
With the pulsational phase at the two observations being so sim- 
ilar, it is hard to imagine that the photospheric temperature could 
be different by 200 K, even for a Mira variable. Assuming a 
Li abundance constant on short time-scales, it is very unlikely 
that the line would go undetected in the observations of Garcfa- 



Hernandez et al. (2007 1. The atmospheric structure of a Mira 
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Fig. 4. T e g vs. Mb i diagram for our sample stars. Filled and open 
symbols denote objects with and without Tc, respectively. The 
stars S Gru and R Hor overlap at (M boU log T eS ) = (-5V35, 3.46). 
The C-star TW Hor has been excuded from this plot. The labels 
give the Li abundance derived (see text). The solid lines are parts 
of isochrones of 400Myrs age (corresponding to an AGB mass 
of ~ 3M Q ) and 6.3Gyrs (AGB mass ~ 1.05 M©), respectively, 
with solar metallicity and no dust formation from Marigo et al. 
(2008). The isochrones have been cut off once the model stars 
become C-rich. 



model, and line strengths might vary considerably with the pul- 
sation phase, even from cycle to cycle. It would be interesting to 
obtain another spectrum of RS Vir to exclude a variable Li line 
strength, as has been reported for a number of carbon-stars by 
Boffi nltaLl ( fT993l l. 

A portion of the spectrum of RS Vir around the 670.8 nm Li 
line is displayed in Fig. [3] with the Li line clearly present. 

4. Discussion 

In Fig.|4]we show the location of our stars in a T e g vs. M\, Q \ dia- 
gram using the data given in Table [T] along with two isochrones 
from Marigo et al. (2008). Most stars nicely follow a trend of 
increasing luminosity with cooler temperature (higher J — K). 
Stars with and without Tc are marked by different symbols. In 
agreement with earlier studies (e.g. |Lebzelter & Hron 2003 



|Uttenthaler et al.||2007a| ), we find no Tc-rich stars below Mbd = 
-4.8. This is another confirmation of the lower luminosity limit 
for third dredge up predicted by various stellar evolution mod- 
els. As in earlier studies, we also find here a few stars with- 
out Tc above the luminosity limit. We refer to the discussion 
Lebzelter & Hron ( 2003 ) that suggests that these stars might 



have too low an envelope mass for third dredge up to occur. 

The data points in Fig. [4] are labelled with the derived Li- 
abundances from Table [2] Standard theories predict a Li abun- 
dance on the AGB below log e(Li) = 1.5. Only one star in our 
sample exceeds this value, namely V441 Cyg. We discuss this 
star in more detail in Sect. |4.1| For the other stars, there seems to 
be a tendency to higher Li abundance with higher luminosity or 
with lower temperature. This contradicts the observations of gi- 
ants summarised in Mic haud & Ch arbonneau ( 1991} and[Mallik 
( 1999), as they find a decrease in the average Li abundance to- 
wards lower temperature. In Fig.[5]we show our sample, together 



with the results for the giants presented by Luck & Lambert 
( [1982) ) and |Mallik| ( [T999l l in a T eff vs. log e(Li) plot. Our sample 
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Fig. 5. Li abundance against effective temperature for the stars in 
our sample (open symbols) and the giants from the list of Luck 
|&Lambert| ( fT9l32l filled diamonds) and |Mallik| ( [1999] filled cir- 
cles). Large triangles represent upper limits from the sample of 
Mallik, small triangles those from Luck & Lambert. The average 
Galactic value for the Li abundance is log e(Li) = 3.1 (Michaud 
|& Charbonneaul[l99T) . 



extends the previous dataset towards cooler temperatures. This 
extension reveals that the trend to decreasing Li abundance with 
lower temperature stops and may even be reversed at the cool 
end. 

The decline in the Li abundance as a function of tempera- 
ture is typically attributed to two effects ( |Mallik| [T999| ) : first, the 
surface abundance will be reduced by a deepening of the convec- 
tive envelope as the stars evolve towards and up the giant branch. 
Second, mass loss should erode outer Li-rich layers during the 
advanced stages of stellar evolution. 

At the lowest luminosities we find in our sample only stars 
with a strong depletion of their Li abundance (see Fig.|4|. These 
stars also have the highest temperature, thus link our sample to 
the one of |Malhk"| ( |1999| ). A deep extension of the convective 
envelope results in the destruction of most of the Li. Lines of 
this element again appear in the spectrum at about M\, Q \ = -4.7, 
i.e. at a slightly lower luminosity than the threshold level for Tc 
occurrence. This suggests that production of Li sets in at some 
point. This is supported by Li preferentially being detected in 
sample stars with long pulsation periods (P > 280 d) and redder 
IRAS [12] - [25] colours. Based on the derived luminosities, 
HBB may be safely excluded as the production mechanism for 
all but one star. 

|Uttenthaler et al.| ( |2007b| ) found four AGB stars in the 
Galactic bulge with a detectable Li abundance of log e(Li) be- 
tween 0.8 and 2.0 (at an error of 0.4 dex) and an effective tem- 
perature around 3000 K. These results agree nicely with our find- 
ings. Uttenthaler et al. identify extra-mixing (e.g. Boot hroyd| 
|& Sackmann| [1999) as the most likely explanation for the re- 
appearance of Li on the stellar surface. We speculate that the 
same mechanism is at work in the Li-enhanceqjstars of our sam- 
ple. The lowest luminosity of a Li-enhanced star in the sample 
of Uttenthaler et al. ( 2007b i is in good agreement with the lower 
luminosity limit found among the targets studied here. 



3 We use here the term 'enhanced' to express a Li abundance above 
our detection threshold corresponding to log e(Li) = 0. This has to be 
distinguished from Li rich objects with an abundance exceeding the av- 
erage Galactic value. 
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In this sample we cannot distinguish whether this change in 
Li abundance is coming from a difference in mass or as a result 
of stellar evolution or both. However, when combining our find- 
ings with the ones from [Uttenthaler et al.| ( |2007b| l, there seems 
to be a lower luminosity limit for any Li production mechanism 
to set in. The existence of such a limit would favour an internal 
mechanism over an external one (e.g. surface contamination by a 
companion - see Uttent haler et al.||2007b| ). We note further that 
mass loss increases substantially during the AGB phase. As mass 
loss should reduce the surface Li abundance, the detection of an 
increase in Li abundance on the upper AGB strongly suggests a 
production of this element inside the star. 

Concerning a possible relation between the Li abundance 
and the occurrence of Tc, we find that all Tc-rich stars show 
also at least a mild enrichment of Li relative to the values found 
for the lowest luminosity stars of our sample. We also find Li- 
enhanced stars without Tc and luminous stars that show no de- 
tectable abundance of neither of the two elements. Furthermore, 
|Uttenthaler et al.| ( |2007b| l also found Tc rich stars with no Li- 
enhancement. This mixture of possibilities indicates the involve- 
ment of a further parameter besides the stellar luminosity. As we 
have no indications of a strong scatter in metallicity in our sam- 
ple, we may assume that mass is the discriminating factor. There 
are some indications that extra-mixing occurs more effectively 
among stars of lower mass, which may lead to an enhancement 
of surface Li in these objects (e.g. Cantiello & Langer 20081. 



If we want to compare our results with the work of Vanture 
et al. ( |2007 1, who find objects with all four combinations of Tc 



rich, Tc-poor, Li-rich, and Li-poor, we face the problem that 
Vanture et al. define stars as Li-rich if log e(Li) > 1.0. All other 
stars are defined as Li-poor. With this approach, only two stars, 
RS Vir and V441 Cyg would be classified as Li rich, with the 
first one at the lower limit. In that sense we would find one Tc- 
rich/Li-rich star (see below), one Tc -poor/Li-rich star, four (five 
if we include TW Hor) Tc -rich/Li-poor stars and a couple of Tc- 



poor/Li-poor stars. This means that all four classes from|Vanture 



et al. (2007 1 would be present in our sample. However, taking 



into account the uncertainties in the Li abundance and the obvi- 
ous scatter of Li abundances between log e(Li) = and 1 , we 
think that the Vanture et al. classes oversimplify the picture. Our 
understanding would certainly profit from finding an abundance 
limit between mildly Li-enhanced low-mass AGB stars and lu- 
minous Li-rich AGB stars producing this element by HBB in 
order to disentangle two different production mechanisms. 

That all four combinations of Tc-rich, Tc-poor, Li-rich, and 
Li-poor are found can be understood at least qualitatively from 
the theoretical point of view. Assuming that 3 He is left from the 
RGB phase, fhermohaline mixing (Cantiel lo & Langer] 2008) 
can connect the H-burning shell of a TP-AGB star with the con- 
vective envelope and thus transport 7 Be to cool layers where it 
may decay to Li that can be observed on the surface. This mech- 
anism is found to be more efficient in low-mass stars (^ 1 .5 M0, 
Cantiello & Langer, 2008| l. This scenario is not constrained to 



stars undergoing 3DUP events, although the faster mixing dur- 
ing 3DUP would ensure that even more 7 Be survives. In higher 
mass stars (> 3 Mq), thermohaline mixing cannot connect the 
H-burning shell with the convective envelope. Nevertheless, a 
7 Be rich pocket can build up above the H-burning shell that can 
be dredged to the surface after a TP, even before formal 3DUP 
occurs. The 7 Be-rich pocket would be dredged to the surface in 
any 3DUP event, which may account for our finding of Li in 
every Tc-rich star in the current sample. The requirement that 
some 3 He is preserved from the RGB phase also introduces a 
mass threshold above which extra mixing can work on the AGB, 



because for stars of too low an initial mass, all 3 He would be 
consumed by efficient extra-mixing already on the RGB. Thus, 
thermohaline mixing critically depends on the mixing history of 
the star. 

In the magnetic buoyancy scenario (Guandalini et al. |2009[ l, 
the mixing is thought to arise from magnetic bubbles with locally 
increased magnetic pressure, thus a decreased density. The mag- 
netic field might be sustained by an a - Q. dynamo produced by a 
region of rotational shear below the convective envelope. Fields 
of different strengths descending from different initial rotational 
velocities may explain the observed spread in Li abundances. 
Independent from 3DUP, Li may be produced on the AGB by 
fast mixing rates (provided a 3 He reservoir is still present from 
the RGB), or destroyed by slow mixing rates. 

Quantitative predictions of models of extra mixing are still 
rather uncertain and require parameter fine-tuning by observa- 
tions. Depending on the adopted efficiency parameter of ther- 
mohaline mixing, Cantiello & Langer ( 2008 ) predict a surface 
Li enrichment of ~ 0.05 dex during one inter-pulse period of 
a 3 M0 TP-AGB model star. After a number of 3DUP events, 
the surface Li abundance may thus increase above the detec- 
tion threshold in a cool atmosphere. Thermohaline mixing is 
expected to be more efficient in lower mass stars, so the sur- 
face Li abundance could be increased by larger amounts than 
for the described case of a 3 M0 star. It should be noted that 
the interaction of thermohaline mixing with magneto-rotational 
instabilities is important, but only full magneto-hydrodynamical 
simulations can lead to a better understanding of these interac- 
tions. 

The magnetic buoyancy picture, on the other hand, predicts 
an upper limit on Li production at the level of log e(Li) = 2-2.5. 
This upper limit is set by assuming a short phase of very fast 
mixing (a few km s~') of magnetised bubbles, in which no burn- 
ing of the material occurs, and the Li enrichment of the envelope 
is set by the amount of "normal" 7 Be that is produced near the H- 
burning shell. The Li abundances below this upper limit can be 
easily reproduced by assuming some Li destruction by a slower 
transport of the material, such that some nuclear burning occurs 
during the transfer to the envelope. If magnetic fields are the ori- 
gin of the extra mixing, both fast and slow regimes of circulation 
can explain very different outcomes with only a single physical 
mechanism. 

To conclude, super-Li rich stars such as V441 Cyg cannot be 
easily explained by the magnetic buoyancy scenario, but the rest 
of our sample would be described well. To decide whether ther- 
mohaline mixing could be the dominant mechanism, more infor- 
mation on the masses and mixing histories of the stars would be 
required. 

We end this discussion with a closer look onto two of our 
sample stars. As far as we know, the reclassification of W Eri as 
Tc-poor has no great consequences on understanding the nature 
of this star. However, the re-classification of V441 Cyg and R 
Hya as Tc-rich and Tc-poor, respectively, has important conse- 
quences, and we thus devote the following two sections to these 
two objects. 



4.1. V441 Cygni 

Even though the SNR of the spectrum of V441 Cyg is low, it is 
clear from Fig. [2] that this star indeed has Tc in its atmosphere. 
It is thus an intrinsic S-star. As further evidence, its spectrum 
around the three strongest Tc lines is shown in Fig. [6] Also 
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Fig. 6. FOCES spectrum of V441 Cyg around the three strongest 
Tel lines. The spectrum has been smoothed with a 5 pixel boxcar. 
The dashed line marks the laboratory wavelength of the respec- 
tive line. Clearly, all three Tel lines are present. 



Wan g & Chen| ( |2002| l list V441 Cyg as "not a good candidate 
Tc-deficient S star" based on its near- and mid-IR photometry. 

|Vanture et al.| ([2007) investigate a sample of 59 Galactic S- 
stars for the correlation of Tc and Li and find all four possible 
combinations of Tc-rich and -poor, and Li-rich and -poor in their 
sample. However, the group of Tc -poor, Li-rich S-stars is formed 
by only one object, namely V441 Cyg. Van ture et al.| ( |2007] l thus 
call it a "unique" object. There was good reason for doubt be- 
cause V441 Cyg actually has never been checked for its Tc con- 
tent in the literature, hence our FOCES observations. As shown 
above, there is no Li-rich (not Li-enhanced\) and Tc-poor S-type 
star known at the moment. 



We thus show that V441 Cyg is Tc- and Li -rich and that it 
needs to be moved to "group 2" as defined by Vanture et al. 
( |2007| l. They offer two possibilities for interpreting this group 
of stars. The first interpretation is that intrinsic (Tc-rich) S-stars 
with high Li-abundance are intermediate-mass (M > 4Mq) TP- 
AGB stars experiencing HBB, thus producing Li in their atmo- 
spheres and destroying C via CN cycling. This picture is sup- 
ported by the low average Galactic latitude of the stars in group 2 
(14°, excluding two relatively nearby stars). However, the bolo- 
metric magnitude, which is available for only four stars of the 
Vanture et al. (2007 ) sample, instead argues against this interpre- 
tation, since none of them has M\, \ < -6™0. Vant ure et al.| (2007) 
thus suggest that their sample consists of a mix of low- and 
intermediate-mass stars. The second interpretation thus concerns 
low-mass S-stars, which might be the result of extra-mixing pro- 
cesses on the TP-AGB (see discussion above). The existence of 
Li-rich C-stars ( |Boffin et al.| |1993| l suggests that low-mass Li- 
rich S-stars could also exist, and have not yet dredged-up enough 
carbon to reach C/O > 1 . 

Which of the two offered interpretations holds for V441 

Cyg? It is located at a very low Galactic latitude of b 1°, and 

our estimate of its bolometric magnitude (M\, a \ = -5'. n 66) is not 
far below the approximate limit for HBB. Concerning the uncer- 
tainties in determining the luminosity we note that only if the star 
were be a fundamental mode pulsator, we would expect a lower 
luminosity. This definitely argues in favour of the intermediate- 
mass HBB interpretation. However, Garcfa-Hernandez et al 



( |2007) > find that Li-rich Galactic OH/IR stars, thought to be 
intermediate-mass AGB stars with a very high mass-loss rate, do 
not exhibit significant Zr enhancements. On the other hand, the 
s-element Rb is observed to be enriched in these stars ( |Garcia-| 
|2006| l. This might stem from the activation of 
Mg neutron source rather than the 13 C(a,n) 16 



Hernandez et al 
the zz Ne(a,n) 

neutron source in these more massive stars. Thus, enrichment 
in Rb is another indicator of high mass. Unfortunately, we did 
not succeed to derive a reliable Rb abundance from our spec- 
trum of V441 Cyg, because the Rb resonance line at 780.0 nm 
is completely dominated by a blue-shifted circumstellar compo- 
nent that makes it impossible to determine the photospheric Rb 
abundance. Nevertheless, we can interpret this as no enhance- 
ment of Rb in the photosphere of V441 Cyg, suggesting that 
it is too low in mass to activate the 22 Ne neutron source. The 
light curve of V441 Cyg changes with a period of 288 d in a 
semi-regular way, and the IRAS [12] - [25] colour, defined as 
-2.5 log(Fi2/F25), has a value of -1.28. Only three of our sam- 
ple stars are bluer than that, which suggests only low obscuration 
by the circumstellar envelope. The variability characteristics, in 
particular the low amplitude of pulsation at this high luminosity, 
could be interpreted as a consequence of a high stellar mass, but 
is difficult to quantify. 

To summarise, the evidence on the mass of V441 Cyg is not 
entirely conclusive, and we cannot clearly decide on which of 
the interpretations offered by Vanture et al.| ( f2007| l holds for this 
star. There are a number of other peculiar intrinsic S-stars with 
high Li abundance, e.g. R And and T Sgr. The Rb abundance 
has been measured in the former by Wallerstein ( 1992| l, and a 
neutron density at the s-process site of 10 s cm -1 was derived, in- 
dicative for the I3 C neutron source operating in low-mass stars. 
An analysis and comparison of these stars would be be very re- 
warding, which we plan to present in a forthcoming paper. 

Finally, Fig.|7]reproduces a small piece of spectrum of V441 
Cyg around the 670.8 nm Li resonance line, along with the 
best-fitting synthetic spectrum based on a MARCS model at- 
mosphere with T e ff = 3400 K and solar metallicity, assuming 
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Fig. 7. FOCES spectrum of V441 Cyg around the 670.8 nm Li 
resonance line (black graph with squares), and a best fitting syn- 
thetic spectrum (grey graph) based on a MARCS model atmo- 
sphere with T e ff = 3400 K and solar metallicity. A Li abundance 
of log e(Li) a +4.6 has been assumed in the spectral synthesis. 
We did not try to fit the blue wing of the Li line, because we in- 
terpret it to be caused by a Li-rich stellar wind from V441 Cyg. 



log e(Li) a +4.6. There is a prominent blue wing in the Li line, 
which we did not try to fit. This feature may be interpreted as in- 
dication of a Li-rich stellar wind that cannot be reproduced with 
our hydrostatic models. Such a wind feature is somewhat sur- 
prising, given the rather blue IRAS colours of this star. Because 
the temperature is probably lower in the wind than in the pho- 
tosphere, more Li atoms are present in the ground state, and the 
wing is deeper than the photospheric line component. Thus, we 
possibly see V441 Cyg in a stage where it enriches the ISM with 
freshly produced Li. 



4.2. R Hya 

R Hya has already been classified as Tc-rich a long time 
ago, albeit with only a low Tc line intensity (Merrill 1952b[ l. 
Spectrogram tracings of R Hya taken at different phases of its 
pulsation cycle can be found in Merrill ( 1952a| l. These spectro- 
grams cover the Tc I lines at 398.497 and 404.91 1 nm, but both 
are absent at all observed pulsational phases. Also in our high- 
quality, high-resolution UVES spectra of R Hya, we are unable 
to detect any of the Tc lines. We find no increase in ZrO band 
strength, wether in R Hya or in any other Tc-poor star in our 
sample. We thus conclude that, with the spectrograms available 
at the times of Merrill, a precise decision about the presence of 
weak absorption lines was not possible, and a misidentification 
occurred in this particular case that has continued in the litera- 
ture ever since. We find R Hya to be the most luminous Tc-poor 
star in our sample, with approximately 10 000 solar luminosities. 
Because of its Mira type variability with a rather long pulsation 
period, as well as its slightly redder IRAS colours, it is proba- 
bly in a more advanced AGB stage than other suspected truly 
low-mass AGB stars, such as ER Vir and EV Vir. 

R Hya is of particular interest because its pulsation period 



has strongly decreased in time ( Templeton et al 



2005 ). It is thus 



Wood & Zarrol 



speculated that it has recently undergone a TP I 
1981). If this interpretation of the period change is correct, then 
either the most recent thermal pulse was not followed by a third 
dredge-up event, or the dredged-up inter-shell material has not 



yet reached the surface of the star. The absence of Tc in the at- 
mosphere of R Hya, however, does not exclude the possibility 
that a recent TP is indeed the cause of the observed period de- 
crease. 

From detailed AGB evolution models (Crist allo et al.| [2009), 
it is known that a 3DUP episode, if it takes place, follows a TP 
with some time delay. For instance, for a 2 Mq model of solar 
metallicity, the first 3DUP episode takes place ~ 670 years after 
the preceeding TP (S. Cristallo, private communication). Thus, 
abundance changes of s-process elements and 12 C on the stel- 
lar surface also become visible only with a time delay. |Wood 



& Zarro 



1981| l suggest that a TP could have occurred in R Hya 
about 550 years ago. If the decrease in pulsation period observed 
in R Hya was indeed caused by a recent TP, we may be able to 
detect surface abundance changes (e.g. Tc) in the future, savin 
a few 100 years, if this past TP is followed by a 3DUP eventrl 

(2002) suggest that the period 



Alternatively, |Zijlstra et al. 



decrease observed in R Hya could just as well be caused by a 
nonlinear instability leading to an internal relaxation of the stel- 
lar structure. While not detecting Tc does not strongly favour 
one or the other scenario, it indicates that possible TPs are not 
strong enough to drive 3DUP in R Hya. 



5. Summary and conclusions 

We analysed high-resolution optical spectra of a sample of 17 
M-type (oxygen-rich) and one C-type (carbon-rich) AGB vari- 
ables for the presence of Tc and the abundance of Li. We em- 
ployed a flux-ratio method to decide on the presence of the 3 rd 
dredge-up indicator Tc in the stellar atmosphere, and measured 
the abundance of the temperature-sensitive element Li with spec- 
tral synthesis techniques of the 670.8 nm Li resonance line using 
MARCS model atmospheres. Bolometric magnitudes were es- 
tablished from near-IR photometry and pulsation periods. The 
lower luminosity limit for 3DUP to occur, established in previ- 
ous works, has been confirmed. We re-classified three sample 
stars with respect to their Tc content, namely that R Hya and W 
Eri are Tc-poor, and V441 Cyg is Tc-rich. We thus suggest that 
R Hya has not yet experienced a dredge-up of freshly produced 
s-elements, but, regarding the possibility that it recently under- 
went a TP, might have a 3DUP episode in the future. V441 Cyg, 
on the other hand, is an intrinsic S-star. 

Lithium is detected in ten sample stars (with a log e(Li) > 0). 
After combining these abundances with bolometric magnitude 
estimates, we find that Li is more abundant in more luminous 
and cool stars. All Tc-rich sample stars are found to contain also 
some Li in their atmosphere. We thus suggest that this fragile 
element could be replenished in the stellar atmosphere by some 
(extra-mixing) process that is only active in the more luminous 
objects. This could also be connected to a mass effect within the 
sample. The S-star V441 Cyg is found to be super-Li rich, with 
a higher Li abundance than the solar photospheric value proba- 
bly by a factor 1000. The gathered evidence does not allow us 
to clearly conclude whether this object is an intermediate-mass 
AGB star (M > 4 Mq) producing Li by the HBB mechanism, or 
a low-mass star with extremely efficient extra-mixing processes. 
The data obtained in this study may be useful for constraining 
theoretical models of such extra-mixing processes. 



4 For the s-process elements to be enhanced, it takes at least two 
3DUP episodes, one to build up a 13 C pocket in the inter-shell region as 
neutron source, and another one to actually mix the s-process products 
to the surface. C enhancement should be present after only one 3DUP 
episode. 
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